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ABSTRACT 


This thesis describes an alternating renewal process and 
a computer program that simulates a transient alternating 
renewal process. The purpose of the simulation is to test 
analytical approximations for the efficiency of a system 


described by such a process. 
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IT. DESCRIPTION OF AN ALTERNATING RENEWAL PROCESS 


This thesis describes an alternating renewal process and 
gives a computer program that simulates a transient alter- 
nating renewal process. Consider a system that can be in 
one of two states, up or down. The system starts up and 
alternates between up and down, remaining in each state a 
Bendom time. The successive times in each state are jointly 
independent. The sequence of up times are identically 
Sis tulouLed With Gist ri butron ae and the sequence of down 
times are identically distributed with distribution 
F xX. and Ys are the corresponding ee duration of time 
the system is in the up state and down state. The means of 
the two distributions are assumed to be finite and are 


denoted by Wy and Wg respectively. The up rate and down 


rate are denoted by ie and Age and are related to the 


means as 
ge See 
ae Sc¥ol jerseys: 
. (1) 
us = i. = mean time down 
d a 





Figure 1. An Alternating Renewal Process. 








Two quantities of interest concerning anvaltermnatang 
renewal process are the availability and the efficiency. 
The availability, A(t), is the probahility that the system 
willbe up at. time t2) The efficiency, E(t) >) 1S thoseratre 
of the expected amount of up time during [0,t] to t. The 
relationship between E(t) and A(t), which will be derived 
later, is 

ic 


fewer ax 


1B es a | (2) 


An alternating renewal process has two imbedded renewal 
processes. Reference 1 gives a detailed description of a 
renewal process. One process, call it a type 1 process, has 
as its renewals the events when the system goes from the up 
to the dcwn state. This 1S a general renewal process with 
the first interevent distribution F, = F and the remaining 


ih 


interevent distributions all being F 


a O Fas where '‘o' 
denotes convultion. The duration of the qige interevent time 
is Y._, + X,; for Mee ee Se CONC INDedde es renewclrmrEOCeas 
has as its renewals the events when the system goes up. The 
common interevent distribution is F = F O Fy: This will 


be denoted as a type 2 process. 

The expected number of renewals during [0,t] is denoted 
by M(t) and is called the renewal function. If Fis a 
function of bounded variation, that is, the difference of 
two monotone functions, and vanishes to the left of zero, 


then 








Down 





Figure 2. Type 1 and Type 2 Renewals 
and the Interevent Times. 


oO 


F” (s) = f exp(-sx) dF(x) 


is defined as the Laplace-Stieljes Transform of F. If the 
times between renewals are distributed as FlsFreees then 


is 


* 
LS) eae ney cs (3) 


For a derivation of this expression see Ref. l. 
For the type 1 imbedded renewal process with renewal 
Pin CctLON M, (t), 
F’ (s) 
ie 1 - a F* (s) 
u a 
For the type 2 imbedded renewal process with renewal 


function M,(t), 


F"(s) F*(s) 
M3 (s) = —4——__$—___ 
1- F fs) F(s) 


Define the random variable S(t) such that 
1 if the system is up at t 


S(t) = 
2 if the system is down at t. 








Down 





2 
Figure 3. Alternating Renewal Process with Two 
Imbeddecd Renewal Processes. 


Let Ny = number of renewals of type 1, 1=41,2, that occur 


moe (0,t)]. Then 
2 


x el 
S(t) = N, Ny 


foelae 
If the expectation is taken, then 


De ih 
ES Ce) E(N,] - E(N,] + J 


(4) 


M. (t) = M, (t) + 1 = A(t) 


The Laplace-Stieljes Transform of (4) gives 


A*(s) = Mo (s) - M) (s) ial 
* * * 
220) 2S i 
* * * * 
1 - Fs) F(s) 1 - F fs) Fs) 
. a(S) 
Lio = ees (5) 


1 - F’ (s) F,(s) 


Let the random variable Q(t) be the amount of up time 
during [0,t] and R(t) be the expected amount of up time over 


the interval [0,t]. Then 





is 
Q(t) =f S(x) dx 
0 


If the expectation is taken, then 


t 
R(t) EG) 1 — ets. S (x) 
0 


t 
Peet) dx 
0 


from which (2) immediately follows | 
1, 


ji NICO ee 


mie) = —— = Ge) % 


If (5) can be inverted, then availability and, hence 
efficiency, can be obtained. For example, if the distribu- 


tions of the up and down times are exponential with cor- 








responding rates a and dae then 
F(t) =l1- exp (-A_t) 
F(t) = 1 - exp(-dA.,t) 
¢ : (6) 
x a 
F; ts) ~ }, + 8 1 = ud 
i 
gives 
Aas 
- (] 
A..+s 
aie) u 
r r 
(| Agts 
_ a + s 
aie + dg + s 





ae 





Inverting yields 


ha ( ds 
A(t) = —— + |-———— exp [- (A_+A ,) t] 
dgtry Mtr, bh rel 
Use of (2) gives 


t A Aa 
if ENC ae { E Jlewwr-a.g))} dx 
u ld el 
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E(t) 


SG Reg 


lied ECA +A 4) 


A*(s) is not always easy to invert. For instance, if 


the down distribution is k-erlang, then 





k 
r 
* 7 d 
Fafs) i eo 
Bo 
+ 1 - Ats 
18) = pie en 
1 - {2} Jc 
3] 
A Qgts) 


(A +s) Agts)* - Alara‘ 


The inversion of (7) can be obtained in principle but it is 
not easily expressed in closed form. If the down distribu- 
tion 1s log normal, a reasonable assumption for many 
systems, then A(t) cannot be obtained in closed form. 

If the system has been operating for a long period of 
time it will settle down into a steady state. For large t 


an asymptotic approximation for E(t) can be found using 











If Neral alse 
lim 0 eee 
t > + oe S =O, = : = 
and 
* 
af (s) =—S=t = Fo (eles am ado mttone 
= s=0 =o 0 
Letting n = 1, then 
t 
if 1G) hx | = ig 
ert © _ 1 dim 29 | etn 
t > © 1: s > 0+ * 
1 oe es 
Ht Fs) Fy (s) 


Wee of WH Hospital's rule for limits im the indeterminate 


form gives 





c 
Ne rey ax ee 
lim 0 a ar =. Byes 
Me ae Ge ~ g§ > 0+ * * * k 
t = ' ! 
[PF '(s)Fo(s) +P (s)Fa'(s)] 
u 
= 5 a (8) 
Wc 


The derivation of (8) was independent of the type of distri- 
butions. If the system is in a transient state, and equi- 
librium has not been reached, then the asymptotic 
approximation is not valid and other approximations must be 
made. Hence, a need exists for testing the approximations. 
The program described in the next section simulates an 


alternating renewal process for this purpose. 
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EE. A SIMULATION OF AN ALTERING RENNER eC ROCEoo 


This section describes a simulation of an alternating 


renewal process. Three distributions are used. 


Exponential 
ite) = dh Seopa) pees 
= 0 otherwise. 
k-erlang 
ko Ae)” exp (-8 &) 
EG) ee 0 t > 0 
i=0 = 
= 0 otherwise 
Log normal 
i 1 2 
Pato exp [- ma) (Qn x-a)'J] ax t > 0 
0 oxv2T 20 
= 0 otherwise. 


Two combinations of the distributions are used. In one 
Betis ination the up times are distributed exponentially and 
the down times are distributed k-erlang. In the other sit- 
uation, the up distribution is the same but the down 
aestrLibution is tog normal - 

R(t), the expected amount of up time, and E(t), the 
efficiency of the system over the interval [0,t], are the 
output parameters and are provided for ten time intervals. 


Analytical approximations for E(t) are given for comparison. 


ikak 








The input parameters vary with the choice of down 
distribution. There are some inputs common to both situa- 
tions: At the unit time interval, Se the up time rate, 


A. the down time rate, the number of simulations and runs 


d 
desired. If the k-erlang distribution is selected the 
parameter k must be provided. If the log normal distribu- 
tion is used the variance must be given. The computer 
program gives more information on the inputs. 

The program functions generally as outlined in the flow- 
chart. Each simulation has these basic steps. First the 
up cycle is entered. While in the up cycle and upon 
exiting, checks are made to determine if a time period has 
been spanned. If not, the next step is the down cycle. 
While in the down cycle and upon exiting, tests are made 
again to determine whether or not the end of a time period 
has been reached. If not, the program returns to the up 
cycle. If at any check the time interval over which the 
Simulation is to occur has been spanned, and, if the desired 
number of simulations has not been completed, then another 
Simulation begins. After the required number of simulations 
have been completed, the output parameters are computed and 
printed. 
The program was written in FORTRAN IV and used on an 


IBM 360/67 computer. 
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